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ABSTRACT 

It is often .considered that taking into account the temperature, 
hence the compressibility, of the electron gas in a plasma is an 
improvement over the cold plasma model. This however leads to a 
number of peculiar and paradoxical results that show the need for 
some ^caution in applying this model. One result is that the boundary 
conditions for the warm plasma do not reduce to those for the cold 
plasma when the temperature approaches zero* Another is that 
evaluation of the impedance or the radiation from an antenna leads 
to widely different results according to the exact size of the antenna. 
This has led some authors to draw completely opposite conclusions as 
to the importance of acoustic waves . Both of these occurrences can 
be traced to tbe fact that the warn plasma equations are of higher order 
than those of the .cold plasma and that the < extra terms contain a small 
factor of the order of a/c, where a i$ the speed of sound; and cthat 
of light . 

This suggests that the techniques of the singular perturbation 
theory can be applied to these problems.. Typically the cold plasma 
model can be applied over wide ranges of the parameters (position, 
frequencies, angle of incidence) and only over narrow ranges forming 
so-called "boundary layers" does ope need to use the warm plasma model. 
Then again simplified equations, can be used and the solutions matched 
on both sides of the layer's boundary. A number of simple examples 
will illustrate this point of view. The analysis confirms that some 
results ape highly sensitive to the values of some, parameters: wire 

radius or gap size for an antenna* temperature of the medium, and 
incident angle of a plane wave . As a result; -the- corresponding 
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"resonances" cannot be^jqbser^ed inf Ip MTctic^^since *^tha. parameters will 

- 

never be realized exactly enough . The boundary layer can then be 
neglected. 

The description of the temperature effect by using the more exact 
kineti? theory is also discussed. The application of this theory is 
much more complicated, and in, a few simple cases where it has been 
worked out, no significantly new, result has been obtained. 
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1. INTRODUCTION 

Most studies of electromagnetic wave and antenna problems in 
plasmas are based on the cold plasma model which assumes that the 
electrons. in the plasma have no thermal energy and form an incompressible 
fluid. In recent years, there have been many efforts to include the 
temperature effect of the plasma by using warm plasma models, This 
can be achieved in. the frame of either . the fluid description or the 
more exact kinetic theory. The warm plasma model. based on kinetic theory 
is quite complicated mathematically and, therefore, has only been used 
in a few simple problems. 

Because of the inclusion of the temperature effect, it is often 
believed that the warxn plasma model , is an improvement over the cold 
plasma model. However, this has led to a number of peculiar and 
paradoxical results that show the need for some caution in using 
this model. .One result is that the field in the neighborhood of a 
rigid boundary for a low-temperature warm plasma is not a small 
perturbation of that for cold plasma, and consequently dqes not reduce 
to the cold plasma field when the temperature is reduced to zero. 

Another is the evaluation of the impedance or the radiation from an 
antenna. Widely different results are obtained depending on the 
exact size of . the antenna. This has led authors to draw completely 
opposite conclusions as to the importance of the acoustic waves. 

Both of these occurrences can be traced. to the following two facts: 

(i) The (fluid) warm plasma equations are of higher order 
than those of cold plasma, and the extra term contains a factor 

6 = - ( 1 . 1 ) 

c 
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where, a is the speed of sound and c that of light (in vacuum) . Assuming 
the electrons in the plasma to form a perfect gas with 3 degrees. of 
freedoms we have 

a = / 3KT /m (1.2a) 

or . _ 

6 « 2.24 x 10 -5 /T (1.2b) 

where K is the Boltzman constant* m is the electron mass, and T is the 
absolute temperature of the electrons in degrees J (Kelvin) . For 
plasmas encountered in the ionosphere or in laboratories for microwave 
experiments, T rarely exceeds several thousand degrees. Therefore, 
the factor 6 is a very small number. 

(ii) The rigid boundary condition for the velocity V at the 
boundary with normal ft, i.e. , 

ft • V - 0 (1.3) 

is enforced in warn plasma but relaxed in cold plasma. In the fluid 
description of plasma, the governing equations are the usual Maxwell 
equations plus the equation of motion for the electrons . For an 
isotropic plasma, the latter takes the form: 


eE 


cold plasma: 

^ “ : iiom(l + v/w) ; 

(1.4) 


- V7 ' V eE 


warm plasma: 

V + 6 ' -r- = — 

(1.5) 


k Q (1 + ivto) iwm(l + iv/w) 



where (-e) , m, and v are the. charge, mass, and the collision frequency 
of the electron, respectively, and k~ = co/c. Note that (1.5) is a 
second-order partial differential equation for V, while (1.4) is of 
zero order. In the limit T or 6 +0, the order of (1.5) is reduced 
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and, consequently, the boundary condition in (1.3) can no longer be 
satisfied. 

The above two facts lead to the creation of the so-called 
"boundary layers" in the field solutions and suggest the application of 
the singular perturbation method to these problems.^ - Typically, the 
cold warm model can be applied over wide ranges of parameters (position, 
frequency, angle of incidence, etc. ) and the result so obtained 
differs from the, corresponding solution in., the warm plasma model by terms 
of 0(6), which are negligible for all practical purposes. Only over 
narrow ranges of parameter forming boundary layers does one need the 
warm plasma model. . Then, in the boundary layers, simplified equations 
can also be used and solutions matched on both sides of the layer's 
boundary. ‘ 

In this paper we will consider: a few simple examples to illustrate 
this -point of view. To avoid unnecessary mathematical., complications , 
the plasma is assumed to be. isotropic throughout. In Section 2, the 
singular perturbation method .is applied to a textbook-rtype problem 
for the purpose of demonstrating the formation of a boundary layer for: 
the acoustic wave in the immediate neighborhood , of a rigid boundary. 
Inside the layer, the acoustic wave contributes significantly to the . _ 
total field so that the .latter may satisfy. the boundary condition in 
(1.3). Outside the layer, it suffers heavy attenuation due to the 
collision loss in the plasma and, therefore, can be safely ignored. 

The singular perturbation method uses this fact specifically and, 
therefore, may simplify the mathematical. manipulations in the warm 
plasma model when applied to sophisticated problems. 

In many boundary value problems In plasmas, we often are 
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interested in physical quantities, such as the reflection coefficient, 
power ratio, impedance, etc. In Section 3, we show with two 

examples that these quantities obtained in a (low- temperature) warm 

- < • 

plasma are nearly the same as those in cold plasma except in , 

boundary layers in parameter space . In these boundary layers, the 
warm. plasma solutions vary so rapidly in parameter space that they 
cease to be physically meaningful quantities. . An attempt to improve 
the situation is to use a more sophisticated model to describe the 
temperature effect in the plasma. Thus, in Section 4, we turn to the 
kinetic theory. However, at least in the example we considered 

/ i . 

(namely, the reflection from a, plasma half-space), the same boundary 
layer, exists in both the fluid and the kinetic theory description! 

In addition to the boundary value problems, another class of 
important plasma problems is the radiation from a prescribed current 
source J(r). An interesting question is whether the fields computed 
from cold plasma and warm plasma (fluid or kinetic theory) are 
significantly different (an indication of the importance of acoustic 
waves). This subject will be discussed in Section 5. The conclusion 
is that as long as J(r) is smooth enough this difference is negligibly 
small, otherwise any sensational results are possible. 
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2 . ATTENUATION OF THE ACOUSTIC WAVE 


The problem under consideration is the reflection of a plane 
wave from an infinite conducting plane in an isotropic lossy warm 
plasma (Figure , 1) . Let the incident field be an electromagnetic wave 
with 

( 2 . 1 ) 


(i) -ik /eT cos8 x] 
H = e ° 1 0 

y 


ik /c7 sin9 z - iwt 
o 1 o 


where. k Q ■ co/c, *1 - (Wp/u)^(l + iv/w) \ Im /e^ >0, 
frequency, and v “collision frequency. The problem is to determine 
the scattered fields which satisfy the wave equations. 


= plasma 


9 x 


,, + k 2 e.cos 2 0; I H ■ 0 
2 o.l o I y 


( 2 . 2 ) 


, .2 2 

a 2 k e . t 

+ ° 1 I P 

a 2 A 2 

9 x 6 


(2.3) 


where 


1 + i - I - 6 2 sin 2 0 

03 , O 


6 = - = ratio of sound and light speeds 


(2.4) 

(2.5) 


P “perturbed pressure in plasma 


and the appropriate boundary conditions at x = 0 and x -*■ «. Since the 
scattered fields will have the same z-variation and time convention 
as the incident one, the common; factor as appeared in [ . J in: (2. 1) for all 




Figure 1. Reflection of a plane wave from a conducting 
plane in a lossy warm plasma 



I 

I 

I 


Figure 2. Reflection of a plane wave from a plasma half space 
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field quantities will be dropped hereafter. 

This simple problem, of course, has an exact solution which can 
be obtained easily. However, in order to demonstrate the existence 
of a boundary layer in the scattered field (as 6 '-*■ .0) in a convenient 
manner, we will attack this problem by a singular perturbation technique 
The basic idea of this technique lies in .obtaining two asymptotic 
expansions for the scattered fields for the smallness of 6. One of 
the .expansions is valid in a very small region cl6se to the conducting 
plane, known as the boundary layer; the other is outside of the layer. 

By matching these, two expansions at an overlapped region, we may 
determine certain constants that. appear in , the expansions, and. finally 
derive an expression uniformly valid inside as well as outside the 
boundary layer. In the present, simple problem our uniformly valid 
asymptotic expression is actually the eract solution. 


To apply the singular perturbation method, first let us consider 
the field expansion valid inside the boundary layer -which is associated 
with the limit 


- cv a. 

X » — fixed, 6 +0. (inner limit) (2.6) 


In terms of the new variable, the two wave equations in (2.2) and (2.3) 
become 


j + 6 2 (e T cos 2 6 ) 
9X 2 1 


H' 


(2.7) 



P p 0 


( 2 . 8 ) 
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Consider (2.7) first. We are looking for an asyniptotic solution of 
the form 


H y * f Q (X) + 6^$) + 6 ? f 2 (X) + 0(6 3 ) . (2.9) 


Substituting (2.9) into (2.7) and equating the terms of same order of 
magnitude, this procedure results in 


9 2 f 


^2 

9X 




9 2 f 


- — — + (e.cos 2 6 )f • 0 . 

^ 1 o o 


When the above differential equations are solved and their results 
substituted intp (2.9), we have 

H y * (C 1 + C 2 h + 6(C 3 + C 4 k) - 6 2 (e 1 cos 2 e o )(C 1 + C 2 j-) 

+ 0(<S 3 ) . (2.10) 

The constants, C's, will be determined later when the boundary 
conditions are applied. The wave equation fbr the acoustic wave in 
(2.8) is a regular one with no appearance of the perturbation parameter 
6, and its solution is simply 

i/e7 t -i/e7 t 

P = Be + De . (2.11) 

To determine a part of the constants in (2.10) and (2.11) , we will now 
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apply the boundary conditions at x = 0, namely. 


E (t) 


(x = 0) = 0, 



0) 


0 


( 2 . 12 ) 


where the superscript (t) signifies the total field, the incident 

plus the scattered. In terms of and P, the two conditions in (2,12) 

become 


6(/eT q sin6 P + i/eT cos9 ) 

1 ° 1 0 S* 




(2.13a) 


to 2 


/e- sin0 

1 Q 

1 + i(v/co) 


(! + H y ) 



= 0 

& = 0 


(2.13b) 


where: . q = — - f . — jr-r- 
M tom(l + iv/co) 

e = magnitude of the charge of, 
an electron 

m = mass of an electron 


Substitution of (2.10) and (2.11) into (2.13) leads to, 


a) , 2 . 

t 6 sin6 (— ^) (1 + C n ) = iqr (1 + i -) (B + D) 

O (jl) ± to 

c 4 

1 - iq(B - D)tan6 

i/eT cos 6 ° 

1 o 


(2.14a) 


(2.14b) 


(2.14c) 


which are four conditions for the six undetermined constants. 

Next, we will consider the field expansion valid outside the 
boundary., layer .which is associated with, the limit 


P 
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X = k Q x fixed, 6 0 . (outer limit) (2.15) 


In terms of the new variable X, the wave equations in (2.2) and (2.3) 
take the form 


3 , 2 - 

— 7 ? + e, cos 0 

9X 2 1 0 


/ 


H = 0 

y 


(2.16) 


2 9 2 2 

6 ; — +tT 

ax j 


p = o 


(2.17) 


and the boundary condition is 


H , P + 0 as X + « 

y 


(2.18) 


This solution of (2.16) and (2.18) is 

i/eT c °s 6 X 

„ . 1 o 

H = Ae 

y 

and that of (2.17) and (2.18) is 


, Im /e"^ > 0 


(2.19) 


P VO 


( 2 . 20 ) 


which means that the order of P is smaller than any algebraic power 
of 6 . 

Summarizing the results obtained so far, we have the inner 
expansion in (2.10) and (2.11) and the outer expansion in (2.19) and 
(2.20), There are a total of seven constants in A, B, D, and C's in 
these expansions, and four conditions in (2.14) for their determination. 
The final step in this method for solution involves the matching of 
these two expansions, which means roughly that the .inner expansion as 
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X -*■ 00 and the outer expansion as X‘+ 0 should be in agreement. 

To carry out such a matching, we need to "stretch" the regions 
of validity for our inner and outer expansions so that they become 
overlapping. Introduce a, new variable 

X n = f (fixed) (2.21) 

n n . : ' 

where n is a function of the perturbation parameter 6 , and is 
asymptotically larger than, 6, explicitly 

6, n 0, but (n/6) . (2.22) 


Note that X = (n/<$)X ■* » and X = nX -> 0i Using limits in (2.21) and 

■ n n 


(2.22) in (2.10), we have 


Inner: H 'v C. + n(C,X ) - n 

y 1 4 n 


:\ 


e,cos 0 « 


+ 0(n J ) 


(2.23) 


where , (2.14a) has been used . , Application of the same limit 
in (2.19) leads to 


Outer: H = Ae . 

y 


t\(±/~z7 cos 8 X ) 

.1 .Of). 


^ A + n (i/e7 cos 9 AX ) - n 
1 o n 


2.l e l cos?9 -o ^2 


+ o(n ) . 


(2.24) 


Comparing (2.23) and (2.24) , we have 


Z, = A, C, = i/eT A cos0 
14 1 < 


(2.25) 
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The matching of P in (2.11). and (2.20) in a similar. manner leads to 

D = 0 (2.26) 

The three additional conditions in : (2.25) and (2.26) , together with 
those in (2.14), determine the two expansions completely. The final, forms 
of the two expansions are 


Inner: 


H- 'v.A 

y 


1 + 6 (i/e"^ cos6 o X) + 0(6^) 


i/e7 xX 
P = Be 1 


i/e^ cos6^ X 

Outer: H = Ae 

y 


(2.27) 

(2.28) 


(2.29) 


P ■v 0 (2.30) 

i . 2h cote 

where: A = ■ . • ,B ■ & -2— (2.31) 

1 + 6h iq (1 + 6h) 

2 1 sin6 tarie , 

; o o 

1 + i(v/w) x , 

{ . : 

f\, . . ’ 

The inner expansion is valid when X defined in. (2.6) assumes a fixed 
value, and the outer expansion when,X defined in (2.15) assumes a 
fixed value. From the .results in (2.27) through (2 . 31) , we may 
derive the expressions for other field components. In particular, 
we are interested in the normal component of the velocity, which is 


h = 


OJ 

_£ 

U) 


found to be 


ik q sin0 

o 

Inner: V v A 

X - me/iT 

o 1 


1 + <$(i/e^ cos0 q X) + 0(6)‘ 


- 2e 


i/^xk 


(2.32) 
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Outer: 


„ l/eT cos0 X 
„ „ lk q sa.n0, . 1 o 

V ^ o n o Ae 
x 


(2.33) 


cue / e t 

o 1 


The final results in (2.27) through (2.33) will now be examined. 

(i) Inside the boundary layer . The acoustic wave plays an 
important part in the total field solution, so as to insure the 
satisfaction of -the boundary condition at x =-0. Away from the boundary 
at x = 0, the acoustic part deqays exponentially as . [see (2.32) ] 

±m</T^ t) "I 


exp 


-X Im(/e~ t) 


exp 


-k x 
o 


6 


(2.34) 


Thus, we may define a "skin depth, ” ; which is a measurement of the 
thickness of the boundary layer, namely. 


6 


\ T) . k i c *i + i h 

O ... 1(0 / CO 


+ 0 ( 6 2 ) 


(2.35) 


At x = X q , the acoustic wave is reduced by e ^ or 37 per cent of its 

2 

magnitude at x = 0. If terms of 0(6 ) are dropped, the skin.;., 
depth is independent of the incident angle 0 q , and, in general, is 
a very small number* . As a numerical example, we have computed Jc fpr 
three atypical sets of parameters encountered in the ionosphere . ; They 


are displayed in Tables I and II. In Table II, we note that 

exists a low-frequency and a high-frequency linlit 

a/w , to -»■ 0 
P 


there 


(2.36) 


2a/v 


u) “ 


which are independent of u>. The numerical data reveal that the skin , 
depth for the acoustic wave is very small except for one case discussed 
below. At high frequency in, the,. F layer, the acoustic skin depth is 



TABLE I 


TYPICAL PARAMETERS IN IONOSPHERE 


Layer 

Height (km) 

T 

a(m/s) 

V 

0) 

P 

D(day) 

60 

300° 

1.16.xl0 5 

io 7 

6x10 3 

E (night) 

90 

200° 

9. 5x10 4 

7xl0 5 

VO 

o 

t — 1 

F 

300 

2DO0° 

3xl0 5 

10 3 

6xl0 7 


TABLE II 


SKIN DEPTH OF ACOUSTIC WAVE x (IN METER) 


6xl0 3 


-2 

r i r\ *• 


** 4 

X 1.9x10 

o 


10 6 

O 

6x10 7 

10 8 

5 .4x10" 2 

2. 6x10" 2 

2. 3x10" 2 

2. 3x10" 2 

*1.6xl0 -1 

2. 7xl0 _1 

2. 7xl0 -1 

2. 7x10" - 1 

5x10" 3 

5x10" 3 

*1.7 

4. 8x10 2 

1.9* ** 10 3 

1.9xl0 2 

3. 1x10 2 

1.9x10 


* at plasma frequency, 

** free space wavelength in meter 
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large because of the small collision loss [see (2.36) ] . However, 

2 

as long as ..(<»> /«)' 5 - << 1, the electromagnetic wave and the acoustic 
wave are practically uncoupled and hence the skin depth is no longer 
a meaningful quantity for measuring the importance of the acoustic 
wave . 

(ii) Outside the boundary layer . The contribution .from the 
acoustic wave becomes negligibly small (due to heavy attenuation) , and 
the total field is nearly entirely made of the electromagnetic wave. 
The amplitude of the 'reflected electromagnetic wave for is given . 


by: 


1 - 6h 


o.l + 6h (2.37) 

where h is defined in (2.31), and A q is the amplitude when a cold 
plasma model is used. Unless the angle of incidence. is Pearly 


(ir/2) , A in (2.37) can, be expanded as 




2h6. + 0(6 ) 


(2.38) 


* Note that, for 0^ = (tt/2) - A with A « 1, we have 

1 - (6 / A ) (to /w)^/ / 1 + iv /to 

P — 

1 + (6/A) (to /co)^/ / 1 + iv /to 
P 

If (6/A) assumes a fixed value, the value of A may differ considerably . 
from unity-. Thus, as a function of 6^, A may be regarded as having a 

7T 

boundary layer defined by 0 < (^ - 0 q ) < A in parameter space where the 
warm and cold plasma solutions do .not agree in the low temperature 
limit. This point will be pursued further in Section 3. 



16 


Thus, A is slightly smaller than’A^. The. decrease , in A is due . to the 
conversion to the acoustic waves which is heavily damped-out due to 
collision; however, one must realize that the .amount of decrease is in . the 
order of 6, and is negligibly small. Therefore,, we may conclude that 
outside the boundary layer, the scattered field obtainable.:.from a .warm 
plasma model is practically identical to that from a cold plasma model. 

(iii) Uniformly valid expression . • Our. inner and outer expansions, 
may be combined to give an expression uniformly valid for all x. The 
standard procedure is to add these two expansions and subs tract out , 
their common part. Take. (2. 32) and (2.33) as an example; the common 
part ( cp ) is 

cp = 1 .+' 6(i/e^ COS 0 O %) + 0 ( 6 2 ) . 

* 

Thus , the . uniformly, valid expression -for V is given by 

ik q sin0 | i/eT cos0 X i/eT t& \ 

V = ° A [ e ■ 1 G -2e 1 .(2.39) 

x we /eT \ / 

o 1 

f ! . ' 

• 

Not surprisingly, (2 .39) turns out to be the exact solution for the 
present .problem. We emphasize the fact that in problems where the , 
exact solution is not available, the singular perturbation procedures 
as illustrated above may often be useful. 

(iv) Lossless warm plasma . As may be seen from (2.35) , in the case 
of w > Wp, the formation of a boundary layer lies in the inclusion of 
collision , loss , in. the plasma model, Even though collision should be 
inevitably present in a realistic plasma, many plasma problems are 
investigated by using a lossless model. The absence of a loss mechanism 
in such an idealized model leaves the acoustic wave unattenuated when . 
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<jo > ojp , and consequently it plays an integral part in the total field 
solution everywhere. Take in (2.39) as an example. When the 
distance along x is measured in terms of the usual free space wavelength 
X Q , the total field, solution has an extremely rapids varying part due 
to the unattenuated acoustic wave. Thus, variation of' a fraction of. 


one per cent in (x/X q ) may result in a significant change in V^. •, This, 
of course, does not agree with our observations in plasma experiments. 
This brings out the point that in using a warm plasma model certain , 
field solutions (e.g. V^) may critically depend on the collision loss 
even though the collision is small. In these situations, the lossless 
idealization may not be a justified one. 



3. BOUNDARY LAYER IN 'FI ELD SOLUTIONS. 

The acoustic wave in plasma, as illustrated in the previous 
section, is heavily damped due . to the collision and, therefore, 
cannot travel a great distance away from the- boundary (or source)* 

Then the next question of interest is that, because of the energy loss 
in the acoustic wave, what is the modification on the the field 
solution for optical waves? We will use two examples to illustrate 
the answer that, except for a narrow range of parameters (boundary 
layers), the field solution for the optical wave in warm plasma differs 
from that in cold plasma only by terms of order 6.. Thus, for most plasma 
problems encountered in practice, the difference is negligibly small. 
However, for parameters falling inside the boundary layers, the -field 
solution for the optical wave may be drastically modified and, furthermore, 
it may vary significantly over a small fraction of an optical wavelength. 
For example, the evaluation of the conductance of a delta-source 
excited cylindrical antenna in warm plasma depends on , the magnetic field in 
boundary layer. Consequently, the conductance may --have .significantly 
different values depending on the precise width of the feed gap. 

In. both, of the examples to be presented below, we will not include 
the collision effect in the plasma for two reasons. First, the 
presence of a small loss modifies the boundary layer somewhat but does not 
change its basic structure. Secondly, in analyzing many boundary value 
problems, particularly the one of antenna impedance, it is traditionally 
based on a. lossless model. By not including the collision here, we 
perhaps can illustrate better why. a variety^ of different results on 
antenna impedance can be obtained depending on the precise width of the 


19 


In the first example we considered the reflection of an incident 
H-wave (H normal to the plane of incidence) from a warm plasma half 

2 

space (Figure 2). The reflection coefficient is found to be 


r — r m . 

o 


(3.1) 


Here IV is the reflection coefficient for cold plasma, 


V ^- 2 


sin 0 - ecos0 

o o 


(3.2) 


-jr. 


2 

sin 0 + scos0 

o l o 


and M is the modification term due to the introduction of the acoustic 
* 

wave 


M 


1 + 5h 


1+ 5h, 


(3.3) 


where 


(1 - e)sin 2 0"” 


0?e- J 2 


sin 0 + ecos0 

o - o 


yfz - 6 2 sin 2 6 


e ■ , 1 - (Up /w) 


Now let us examine H as a: function of e and sin 0^. For a given 


6 and 6 -*■ 0, it may be shown tjiat 


M = 1 + C) (6) 


(3,4) 


* To include the collisioii loss, simply make the following two replacements 

(w /a )) 2 .. , . 

e = 1 -_JL_ , 6 - ? - 1 ' . 


(1 + iv/oi) 


/ 1 + iv / c 
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for all 0 - sin 2 0 Q - 1 and e. <1 except when 


sin 2 0 = -r— | — • + A 
o 1 + e 


(3.5) 


where A is a small number. With A =0, (3.5)gives the condition for 
total transmission in cold plasma (F q ■ ,0, but not T) . We will now 
examine the expression for T in (3.1) under the condition of (3.5), 
that isj the boundary layer of the reflection, coefficient. As indicated 
in Figure 3, there are two subdivisions depending on whether e itself 
assumes a fixed value or. a value comparable to 6. 


(a) s is fixed . One may- show 

2 


r ^ 

o 


(1 + e) (s - 1) 
As 2 


M <v'l + (— ) — - — ( — £ — ) 3 / 2 
n 1 + V 1 + S K 1 + e } 


Thus, the reflection .coefficient- for warm plasma is of 0(e) when 
(6/A) -*■ °°, or of 0(A) when (6/A) ->■ 0. In either case, the difference 

between T and T is negligibly small. 

(b) e is comparable to 6 . The situation is more complicated. 

Let us concentrate on the case of . e « 6 ; 5 and A are the variables. 

» i , . . • ] 

The main results of our study are summarized in Table III where the 

values of F o , M, F> R> and n for different orders of (A/6) are listed. 

The parameter n is, the , ratio of acoustic power in the plasma end the 

incident power from the free space. The numerical values in Table HI 

indicate the extremely rapid variation of, the field solutions in 

2 2 _5 

the boundary layer. For example, when sin 0 q “ 6 = 10 (corresponding 

to A = 6^; T = 2000° and 6 = 10 "*) , the reflection coefficient ,F is 
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TABLE III 

FIELD SOLUTIONS IN BOUNDARY LAYER 


i 

A i 

1 

r 

o 

M 

r 

n 

0 ( 6 5 ) 

0 

\ 

00 

1/3 

4/9 < 

6 4 1 

-l 

0 

0 

1 

G(6 3 ) 

l 

1 

1 

0 

) 


zero and the total incident power is totally pumped into the acoustic 
waves (o = 1) . Thus, it can be said that the warm plasma model gives a 
completely different result frbm that of the cold plasma (r Q ■ -1) . 

However, if we vary the incident angle by 5x10 ^radians so that 

2 2 3 ■ 3 

sin 0 o =6 + 6 (corresponding to A » 6 ) , the ref lection coefficient 

T becomes unity, in agreement with the result in cold plasma. Thus, 

total transmission .and total reflection are separated by merely 

5 xlO 7 radians of the incident angle! 

In conclusion, . in the problem of reflection from a warm plasma 
half space, the field solution is : practically the same as that , 
obtainable from a cold plasma model except when the parameters (e, 9 q ) 
fall inside a boundary layer. As indicated in Figure 3 (not to scale), 
the boundary layer is a region around e = 1 - .X, Wp/w) ■ 0, having a 
dimension in the order of 6 in the parameter space (e - 0^ plane) , The 
field solution for parameters in the boundary layer varies extremely 
rapidly, and, therefore, has a dubious physical significance. 

In some electromagnetic boundary value problems, the physical 




22 



sin 2 0 



Figure 3. In the dotted region (not to scale), the cold 
and warm plasma solutions differ significantly 
for the problem sketched in Figure 2 



Figure 4. An infinitely long cylindrical antenna in a 
warm plasma 
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quantity of interest falls exactly within 5 the "boundary layer of the 
parameter space, and this may lead' to a variety r of different conclusions. 
Such situations will be illustrated in the second, example* namely, the 
calculation of the impedance of an infinitely, long cylindrical antenna 


immersed in an isotropic lossless warm plasma.... -The geopetry of the 
antenna is shown in Figure 4; it is excited by- a delta source with 
unit voltage amplitude located at z ■ 0. The. expression for the 


induced current on the antenna has been obtained by a number of 
3 4 

authors * and is duplicated below 



iaz , 
e da 


w 


' il o 1 '( ? e a) ' 

r _2 . + 

6 H^C^a) 


2 2 H (1) (£ a) 

’ a ' -o-, p - 

S H. (1) (? a) 


(3.6) 


p -< 


where 


? = /. 2 2 = +i/7T 72 n ®;;e,p 

n k - a a - k , 

n n * 

k e ■ >0 ’'I - (» M 2 • V" /c .- 

P ' 


k ■ k /5 , 6 = ratio of acoustic speed and 
^ e light speed 

The integration contour in (3.6) is slightly above. the real axis for 
a < 0 and below for a >0. The admittance of the antenna defined 
as 


Y = G + i? .= 1(8 - z .) (3.7) 

- ' _ O 

where z q is a suitably spall! distance from the idealized feed at z = 0 
and may be identified with the half-gap width of the actual feed. We 
recall that even for an antenna situated in the free space, the imaginary 
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part of I(z) for z close , to the gap has a logarithmic singularity, and, 
therefore, has , a boundary layer. In order not to. confuse that boundary 
layer (because of the delta source) with the boundary layer to be 
discussed below (because of the acoustic wave) we will concentrate 
on ,the real part of I(z) or the evaluation of —the conductance G (for 
a) > u )• • Furthermore, to facilitate the estimation of certain integrals, 
we will assume that the parameter 


/ " • 2 
A = k a = .k a.l - (w /w) 
e o p 

is reasonably large so that (A) can.be well-approximated by its first 
* 

asymptotic term. Summarizing, our, task is to evaluate the real part of 
I(z = z q ) given in (3.6) under the assumptions 

6 0 , A. >> .1 (3.8) 

There are two subdivisions, depending on the relative magnitude of 
z^ and the acoustic and optical wavelengths 


„ k z / r fixed (inner limit) (3.9) 

CD 4 - k , - -V 2 'i- - (»>) - 

p O 6 p 

2 

(ii) Z « k g z ® fixed (outer limit) . (3.10) 

We will first concentrate on, the inner limit, where the gap width is 
so small that it has to be measured in thrms of the acoustic .wavelength. 

m 

Under the assumptions in' (3. 8) and (3.9), the real part of I(z q ) in 
(3.6) may be approximately evaluated with the result 

* 

G y G q M (3.11) 

* The error is less thah 5 per cent when A = 3. * ' 
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Here G q is the well-known (approximate) cold plasma solution for a 


thick antenna 


k a 

, _ 

*o 120 


(3.12) 


and M is the modification factor 


M( 2 () = 1+2 (s' - 1 ) COS (a 2 ). + 


,^2 r.fa/Z cps(Zsin 


* 7 0 y (u/w ) 2 + tan 2 6 


d9 


+ 0 ( 6 ) 


(3.13) 


The second term in (3.13) Is duk to a surface wave contribution coming 
from a simple real zero of the denominator of the integrand in (3.6a). 
Under the condition in (3.8) , the surface wave has a propagating constant 
(along z) approximately given by 


a ^ k s 
P 


VS 


(u p /u) 


(3.14) 


4 - 

provided that is not too small . The- third term in (3.13) 

comes from the portion of the integral in (3.6) in the range 

k < a < k (3.15) 

e ' P 

In the cold plasma model, there is -no surface wave an4 there is no 
contribution to G from the integral range k^ < a (invisible range) . Thus, 
the deviation of M(2) from unity is the special feature of the warm 
plasma model. The following observation is made regarding M(2f) in 
(3.13). 

: ‘ j ■ 

(i) M does not reduce to unity as 6 ■+ 0. Thus when the gap width 

z q falls inside the boundary layer [i.e. (3.9)], the cold plasma solution 
G cold cannot recovered from the .warm plasma solution G by letting 
the, temperature (or 6) go to zero . The same conclusion holds even if 



26 


a small loss, in the medium is introduced* 

(ii) M is an extremely rap id- varying function of k z = k. z"\/l-(w /w) 

e o o ov p 

(the gap width measured in terms of the optical wavelength in , the plasma). 

This follows from the fact that ^ * k z /6 . and ld. -» 0. Thus, the 

' e o . • 

conductance of the antenna, depends critically on the precise gap width 
in terms of the optical .wavelength. 

Next, let us consider the outer limit ; in -'(3.10) . It may be shown 
that the expressions in (3.11) and (3.12) are still valid. Now note that 
the third term in (3.13) behaves as 


2 

7T 


fir/2 



cos (6 Zsin8)d9 . 

(u/u + tan^Q . 
P 


0 ( 6 ) 


(3.16) 
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H(Z) =1+0(6) 


for lossy plasma. Thus, an agreement between .the warm and cold plasma 
is obtained. 


The conclusion of this antenna problem may be stated as follows. 
The question of interest is whether the conductance of an infinitely 


long cylindrical antenna in a warm plasma is a perturbation of that in 
a cold plasma, and whether these two solutions agree with each other 
in the zero- temperature limit . Our analysis shows that the answers 


depend on z q , the width of the feed gap. If the gap width is small 
compared with the optical wavelength (so that z q is fixed), .the 
warm plasma solution cannot be reduced to the cold one when. the 
temperature approaches zero. Furthermore, the conductance in warm 


plasma is an extremely rapid-varying function of k^ z q , and, therefore, 

is not a physically well-defined quantity. If the gap is comparable 

to the optical wavelength and if the plasma has small collision loss, 

the Warm and cold solutions for the conductance then become almost the 

same. However, in solving antenna problems, k Q z q is usually assumed 

to be so small that k^ z q is fixed. This explains why a variety of 

different values for the conductance can be obtained, depending on the 

exact value of k z . 

oo . 
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4 . THE WARM PLASMA MODEL BASED ON THE KINETIC THEORY 

In addition to the (fluid) warm plasma model used in the previous 
section, the temperature effect can be also described by the more 
sophisticated kinetic theory at the expense of the mathematical 
simplicity. To date, one of the few electromagnetic wave boundary 
value problems that have been solved by using the kinetic theory is the 
reflection from a plasma half space. In Section 3, we have examined 
the solution of the problem based on the fluid theory and have found 
the existence of the boundary layer in .the parameter space for the 
field solution (e.g. reflection coefficient). In this connection, an 
interesting question is whether the more sophisticated kinetic theory 
can rescue us from such a, singular behavior in the field solution for 
a low-temperature plasma. 

As in the fluid model, a crucial boundary condition at the plasma 
free space interface is one regarding the velocity of the electrons. 
Referring to Figure 2, instead of the rigid condition on, the mean 
velocity of all electrons used in the flnid model 

V x (x = ,0,z) = 0 (4.1) 

the condition of specular reflection is commonly employed in the 
kinetic theory, namely, 

f(v_,v_) = f(-v ,v_) at x .= 0 (4.2) 

x z x z 


* Alternative names are "transport equation model" and "hydrodynamic thoery." 
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where f is the perturbed distribution function for the electrons, 
Recalling that 


V x (x ■- Q,z) - 


V f(v ,v ) 
x x’ z 


x = 0 


dv dv_ ; (4.3) 

A 2 


the enforcement of (4.2) implies (4.1). 

Using the linearized coilisionless Boltzman-rVlasov equation and 

the usual Maxwell equations, the reflection coefficient for an incident 

H-wave (H ,, E , E ) from. a plasma half space subject to the condition 
y x z 


in (4.2) is found to be 


5,6,7 


z cos"e - z 

o o 

z cos 2 e + z 
o o 


(4,4) 


Here Z q = , and Z is the surface impedance of the plasma 

half space defined by 


Z = - 


H • 

y 


X :■ .0+ 


(4.5) 


and is given by 


Z = lim 


u 


0+ ilTe o 


, c 2 k 4 D T (k) 


o) 2 k 2 D^(k) 


ik 

e x dk (4.6) 
x 


where k = k x + k z, and k 
x z * 


The two functions D^(k) and D^(k) 


are plasma dispersion relations for transverse and longitudinal waves, 


respectively . They are given by 

v k > - 1 -o 2 - ^ 


f o (v) 


k • v - to 


,2 

d v 


(4.7) 
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D (k) =1 - a) 

P 


f o (v) 


d 2 v 


(k • v - w)‘ 


(4.8) 


where f (v) is the unperturbed distribution function for the electrons, 

° n — t 

and is assumed to, be isotropic (depending on v =\/v + v , not on v) . 
At this point, it is interesting to mention that if we use the plasma 


dispersion relations obtained in the fluid model, namely 

2 2 
ll) —,u 

e t (W • 1 - p y 2 . 


V k > 


1 - 


ck 

2 
w ' 
P 


2 2 . 2 

w - a k 


(4.9) 


(4.10) 


in (4.6), then the integral can be evaluated explicitly to yield the 
result : 




sin 2 0 o + (1 " £) 0 o 

- sin 2 0 
<5 2 


(4.11) 


Substitution of (4.11) into (4.4) recovers .(3.1) , as expected. 

Return to the solutions obtained by the kinetic theory as given 


in (4.4) through (4.8). Before , trying to evaluate (4.6) for some 
assumed f^(v), let us first point out the general features of the present 
solution which are different from the fluid solution 

(i) In addition to zeros, the functions B^,(k) and D^(k) may 
have branch cuts in the complex k x -plane. The location of the possible 
branch cuts is determined by the condition 

k • v - oj “ 0 . . (4.12) 

Since in the present problem k^ = k Q sin9 Q , this condition becomes 
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co - v k ; sin0 
z o o 


(4.13) 


If we assume the unperturbed electron velocity cuts off beyond a fixed 
number v 


f (v) .» 0, for v > v 
o * o 


(4.14) 


then we may show, that the branch cuts in the complex k^- plane for a given 
co (with Imw > 0) is determined by two equations (Figure 5) 


Re k 


M to 
Re co 


(4.15a) 


\\ ■>- 


Is)' eiA 

C 0 


1 - 


=> a . 
o 


(4.15b) 


In evaluating (4.6) by deforming the contour in , the upper half k^ plane, 
we note that the contribution to Z ,not only includes those .from the 
zeros of D„, and D T but also from a branch-cut integral. The addition of 
the branch- cut integral in the kinetic theory accounts for the existence 
of the van Kampen modes (modes with continuous spectrum) in the plasma 
half space. 

(if) In the limit Imoo -► 0+, the dispersion functions in (4.7) 
and (4.8) may assume complex value. Let us concentrate on (4.8) , which 
may be rewritten as 


o .. 


V k) 


1 - uf / F 0 (V 1 ) : . dv. 

p '. v (kv - iiy; 

o 


(4.16) 


where ( v i> v 2) are the components of . v (parallel, perpendicular to k) , 
and 
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The (positive) imaginary ^art in (4.18) indicates that the longitudinal 

wave with wave number |k| > (oj/v q ) suffers a (Landau) damping even in a 

collisionless plasma. An interesting consequence of Landau damping is that 
2 

when (to/u>p) < 1 (cut-roff conation for plasma) , the surface impedance 
Z in (45) still has a real part and hence the total reflection ( | F ] = 1) 
can never be obtained. Similar comments apply to the transverse wave and 
its dispersion relation. 

Whether the above two differences between the kinetic theory and 

the fluid theory can lead to significantly different field solutions 

• ' • ' ’ ; /’ ii- y 

depends largely on the assumed unperturbed distribution function f Q (v) . 

■ ' "i l . . ‘ ’ 

Commonly, the Maxwellian distribution is, used for f^Ov), namely 
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f 0 (v) 


i ' 2 

r • i • 2a 

e 

2*a / 


(4.19) 


where 


h/. 2 , , 2 

=, Vv + y. 


a = sound speed in the electron gas . 


For such a f^(v), the reflection coefficient T given in (4*4) has been 
approximately evaluated by Wes ton His result is given below 


( I’ o + 0(6), if e is fixed 


f w 

a, \ 


\ A + <?g_ i if e. « 1 


(4.20) 


1 + 5 §+ 


where f is the reflection Coefficient obtainable by using a Cold plaSma 
model and is given explicitly in (3*2) * and 


Sin^Q 


(\/e - s in^ 6 „ + e cos 8 ^ - 6^sin^ tk 


J _ I CV.VOU . I W 1 

o - o » 1 - e 


(4.21) 


This result is practically identical to the Solution obtained by the 
field model in Section 3 i Thus* the boundary layer also exists in 
(4.20) , the reflection Coefficient obtained by using the kinetic theory! 


* Comparing (4*21) and (3*3) we note that g_j. = h + + 0(e) 
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5 o RADIATION IN WARM PLASMA 


The creation of the boundary layers in the field solution of warm 
plasmas is attributed to the fact of small .5 and to the enforcement of. 
the rigid boundary condition. An interesting question in this connection 
is whether the warm plasma model can give any significantly different 
result when the rigid boundary condition is not enforced. An important 
problem belonging to this category is the radiation from a prescribed 
current source J(r) in an unboun4ed warm plasma; it will be studied in 
some detail in this section. 

, 8 9 & 

It has been shown by several authors ’ that, the Fourier transform 

of the electric field E(r) , denoted by E(k) is related to J(k) through; 

the general relation 


E(K) = 


i( ^o 


k 2 I kk« : 
k 2 D T (k) 


2rr 
c kk° . 

oo 2 k 2 D L (k) 


J(k) 


(5.1) 


The above expressions are valid for both. (isotropic) cold and warm 
plasmas (fluid or kinetic theory) . For warm plasmas , the explicit forms 


of D„ and D T are given in (4.7) and (4,8) for the .model based on the 
i L . 

kinetic theory, and in (4.9) and (4.10) for the model based on the fluid 


theory. For the expressions of D^ and D^ in cold plasma, we may simply 
set the sound speed,a = 0 in (4. 9 ) and (4,10). 


The result of E(k) as computed from (5.1) may be classified into two 
types , according to the smoothness of J(r). The first type occurs 
when J(r) is smooth enough so tjiat J(k) is nonrzero only -when - (6k/k ) << 1. 
In such a case E(k) in (5.1) depends ..on the values of D^,(k) and D L (k) only 


The-Fourier -trans f orm-is - F (k) - ® 


: i 

r°\ 

1 J 

f 3 

— oo 


F(r)-e : 


-ik°r 


,3 

dr 
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in the range 


(Sk/k ) « 1 
o 


(5.2) 


It may be shown that, under the condition in (5 , 2) , 


D < k > warm ’ D « cold + 0 


1 k6 I 2 


(5.3) 


where D stands for either or D^. Explicitly- we have 


IQ 


cold plasma: 


D T (k> - 1 - r 


(0 

1 _ -P-- 
2 
CD 


(5.4) 


D L (k) = 1 - 


it ! 2 

0) 


warm plasma (fluid theory) 


D x (k) = 1 - 


/ k 




0) 

1- _JL_ 

1 2 

U): 


(5.5) 


D L (k) = 1 - 


to 

,0) , 


i + 


k6 

k 


k<5 


k 


+ 


warm plasma (kinetic theory) : 


T 


Ik 1 


2 .2 /to 

2 


k6 

k 


Mi 

+ T* f j 


1 + 

\ 1 

' 0) . 

\ i 

-l 



i o 


+ 


(5.6) 


/ 0) 
M- 

2 

i + 

k6\ 2 _6_ 

r k6 

tU) 


Ik 3> 

k 

1 / 



l o/ 

[ o 


+ 


In the formulas in (5.6), 6 = a/c and a ="/ 3KT/m is consistent with 
those in (5.4) and (5.5). 
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A study of (5.4) through (5b 6) reveals that 


E(k) l « E(k) 


warm 




(5.7) 


except when 


iLi 

i-U 


0 ) 


0 ( 6 2 ) 


(5.8) 


^-.0(1/6) . (5.9) 

When condition (5.8) is satisfied, D^(k) for warm plasma may be 

& 

significantly different ^from D T (k) for cold plasma; while the, 
satisfaction of (5.9) implies a remarkable deviation between the D^,(k) in 
the kinetic theory .model and D^,(k) in the other two models. Except for 
these two very special cases, we conclude from (5.7) that there is 
no essential difference between the radiation field produced by a 
given smooth current source whether. a cold or a, warm plasma model is 
used. This conclusion is established on a general term with no 
particular reference on the current source, as long as it is sufficiently 
smooth so that (5.2) is satisfied. . 

A special feature of the .kinetic theory model, as has frequently 
been mentioned in the literature, is the fact that D (k) and D (k) 

* ' L . Li t 

defined in (4.7), (4.8), and (4.9) have an imaginary part even for real 

* 

k in a lossless plasma, accounting for the Landau .damping effect. Take 


W -■ . 

Recall the .problem of reflection from a plasma half space discussed 
in Sections 3 and 4. , This scattering problem is equivalent to that 
of a current sheet radiating at the interface. This explains why 
a boundary layer exists in the field solution when (5.8) is satisfied. 
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D L (k) as an example: its imaginary part is given, for real k, by 


Ip D l 00 


U) 

2 

Mc \3 

- 

U) 


1 kfj ex P 

- f <V w)2 


(5.10) 


A direct consequence of (5. IQ) is that when w < the complex power 
radiated from the source 


t - -J e(t) • J(r)* d 3 r (5.11) 

may still have a real part, indicating dissipation in a lossless plasma. 
However, we must note that when the condition in (5.2) is satisfied, the 
imaginary part of D^(k) in (5.10) is exponentially smallef than its 
real part in (5.6). Thus* with^ a lightly lossy plasma [ i*e. , let 

(w /w) (u /u) / (1 + iv/u), and 6 -> 6/ /l + iv/co with v being the 

P P 5 - : - 

collision frequency ], the effect .of the Landau damping is not as 
important as that of the collision. 

The other type of result which may be obtained from (5.1) is the 
one associated with J(k) which does not cut off .fast enough and may 
assume finite values in the neighborhood Of k ■ k Q /6 . In such a case, 
E(k) computed from the warm plasma model may be significantly different 
from the corresponding E(k) computed from the cold plasma* This may be 
illustrated by the following example* 

A popular method in evaluating antenna impedance in plasma 


It is known that the longitudinal wave, not. the transverse wave, 
suffers Landau damping in A plasma. D«*00 also has an" imaginary part 
because the relativistic effect has. nor been included in the Maxwellian 
distribution in (4. 19) * 
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is the so-called "included e.m.f. method" which entails an assumed 

1 t 

current distribution and then computes the antenna -resistance through 
the formula 


R 


7 *•/' 


' / E(r) - J*(r) d 3 r 


(5.12) 


where I is a normalization factor. For the case of a, linear antenna with 

O . ..... i ‘ ‘ v 

radius a and half length h, J(r) is commonly assumed to, be "triangular, " 


J(r) 


v 2ira 


z 


1 - 


6(p - a) z , 


for z < h 


(5.13) 


The Fourier transform of the , current is 


J(k). ■ I 


Q, 1 


sin^khcos^) 
(khcos|). 


J^ka sin8) z 


(5.14) 


where (k, 0 , 4>) are the spherical components of k, and J q is the zeroth . 
order Bessel function. The rat® of cut, off of J(k) for large k depends 
on the value of h and a. In the following, let us concentrate on two 
special cases?. 

(A) k Q h =,0(1) and k, a ■- 0(6 1 ^ 2 ) (5.15) 

(B) k h = 0(6 1 ^ 2 ) and k a » 0(6 2/3 ) (5.16) 

o f "* o 

The choice of the above parameters is designed to make the results 
physically illuminating as well as mathematically simple. For large k 
such that 
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we have 


J(k) 


K » •; — , fixed 
k 

o 


0(6^^), for case A 


0(6) , for case B 


(5,17) 


(5.18) 


(5.19) 


As we will show next, J(r) in case A. is sufficiently smooth so that there 
is no appreciable difference between the cold and the warm plasma 
solutions for resistance ,R in (5.12) , but J(r) for case B is not, Under 
the condition that 

k eh r— . 

-r—— >> 1 and h >> a, where e ■ VI - (to , /to) * (5.19a) 

o P 

an approximate expression for R is found to be 


R - R M 

:• oi 


(5.19b) 


where R q is the resistance for cold plasma and 


M 


,3if(l - e ) 

— ii - - - i 


i% 2 1 

k ae \ 

O 

1 + " . 3/2 
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Then it is a simple matter to verify that 

M ■ 1 + 0(6^ 2 ) i f 0 r the case A 

which is essentially unity. However, for the other. case we have 
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which is a rapidly Oscillating function., If the radius of the antenna 
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is changed by a small fraction 
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where X q is the free space wavelength, M in -(5.21) may vary from unity 

to an extremely large number [in the order of 0(<5 - ^^)]* In other 

words, the acoustic wave "can contribute from zero to nearly 100 per 

-3 

cent of the total radiation resistance by a variation of >10 X or 
10 ^ X in . radius ! 

p 
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6. CONCLUSION , 

The warm plasma model based on the fluid theory has been used 
extensively in the literature. In the last- several years there seemed 
to have been a trend toward treating by this model every electromagnetic 
boundary value problem previously solved by using the cold plasma model. 

The primary motivation for doing that is to study the effect of the 
acoustic wave which is not accounted for in the cold, plasma description. 
However, this effort does not seem to have -led to a definite conclusion 
on the importance of the acoustic wave. Oftentimes, in closely related 
or even identical problems autkors have drawn completely opposite 
conclusions, depending on thdsir choice of parameters. In the present 
study, we have discussed an explanation of this, based on boundary layer 
theory. The main results of our study concerning the (fluid) warm plasma 
model may be summarized as follows: 

(1) In a low-temperature plasma, the acoustic wave constitutes 
an important part of the total ^ fie Id only in narrow-layered regions close 

i 

to the rigid boundary. Outside the layers, even a small collision frequently 
produces an attenuation of the acoustic wave much larger than that of the 
electromagnetic wave. Therefore, the total field consists almost 
exclusively of electromagnetic waves . 

(2) Due to the excitation of the acoustic waves (dissipated 
in collision loss or not), the solutions of physical quantities which 
are of interest in electromagnetic studies are modified. Examples of 
such physical quantities are the reflection coefficient of an incident 
electromagnetic wave, radiation power from a current source, input 
impedance of an antenna, etc. The modification is such that, except 
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for a. few isolated regions in the parameter space forming the so-called 
"boundary layers," the solutions obtained' by using warm and cold plasma 
descriptions differ only by a negl'igibiya small amount., However, inside 
the, boundary layers the warm plasma solutions vary so rapidly in the 
parameter space that they cease to be -physically meaningful quantities. 

(3) The existence of the'boundary layers accounts for an 
important fact that warm plasma solutions may be extremely sensitive tp 
the exact value of ..the parameters. This - perhaps is the main explanation 
for contradictory conclusions about the importance of .the acoustic waves. 
Thus, in the study of electroipagnetic problems by using the warm plasma 
model, a systematic investigation of the ! solution dependence on parameters 
is essential. This is particularly so trtie'n the solution is .not given in 

a simple analytical form and numerical computations are necessary. 

Results based. on spot calculations may be very misleading. 

(4) Those sensationally different results in the boundary 
layers of the warm plasma solution probably cannot be observed in, 

practice since the parameters will never be realized exactly enough. 

( . 

Furthermore, the warm plaBita/ opodei Itself perhaps is not a good description 
of physical plasma when tHosh .situations arise. 

(5) Another attempt to assess the importance of temperature 
effects is to use the more exact kinetic theory. The application of 
this theory is much more cpihjjlicated and in the flew simple cases, 
where, it has been worked out, it- has not r produced .significantly 
new results. In a problem (i.e. , reflection from a plasma half space) 
for which an .explicit solution can be found, the solution has practically 
the same boundary layers in the parameter space as those found in the 
fluid model. Outside the boundary layers, the solution from the. 
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kinetic theory is again approximately the cold plasma result. 

In addition to the boundary value problems, we have also considered 
the radiation in an unbounded plasma from^a given current, a problem 
in which the rigid boundary condition on the electron velocity is not 
applied. As long as the, current is so smooth that it does not vary 
significantly over an acoustic wavelength-,- the radiation field in cold 
plasma and in warm plasma .(fluid or kinetic theory) are practically 
the same., Furthermore, the often mentioned Landau damping in such cases 
is not as important as the other loss mechanism such as collision. 
However, if the current is not smooth, the warm plasma solutions may 
be significantly different and, in some cases, highly implausible. 
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